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Abstract 

Average persistent current over a set of diffusive metallic rings with fixed 
number of electrons is considered. We study the the case where the phase 
breaking time is much greater than an inverse average inter-level distance. In 
such a situation, many return events for an electron have to be taken into 
account. As a result, one arrives at a non-perturbative problem for a fixed 
by an external magnetic field Cooperon mode. This multi-Cooperon problem 
problem has been considered previously by Altland et al., Europhys. Lett., 2, 
155 (1992) and in several following papers within the framework of supersym- 
metric approach. Such an approach involves very tedious calculations which 
were performed using computer algebra package. Here we solve the problem 
in question with the help of replica trick. It is demonstrated that the replica 
trick in combination with a proper analytical continuation in the replica space 
allows one to obtain the result in much more explicit way. 

PACS numbers:72.10.Bg, 73.20.Dx 
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I. INTRODUCTION 



Magnetic properties of small conductors were extensively studied during the last several 
years (see§@ and references therein). It has been understood that the magnetic moment 
(and the associated persistent current) induced by an external magnetic flux is a very spe- 
cific manifestation of mesoscopic behavior. While originally predicted to appear in clean 
one-dimensional metallic rings!, most of the recent discussion about persistent currents has 
focused on metallic rings containing impurities^. Static magnetic properties of small rings 
and dots were studied by several authorJUHI. An important step in the understanding of 
magnetization of mesoscopic quantum rings was taking into account the difference between 
canonical and grand canonical ensembles00&0. It has been shown that the magnetization 
of isolated rings with fixed number of particles is much larger than of the ensemble of rings 
kept under fixed chemical potential. As a result, the main contribution to the magnetic 
moment has been expressed in terms of the fluctuation of the number of particles at fixed 
chemical potential, ((5N) 2 ). The latter quantity has been analyzed in00 under the condi- 
tion h/Ar^ ^> 1. Here A is an average inter-level distance at the Fermi level, A -1 = uV {y 
is the density of states at the Fermi level, while V is the volume), is the phase-breaking 
time. 

Let us discuss the physical meaning of the parameter H/Ar,. As well known§, in the 
absence of external magnetic field the quantum correction to the conductivity is proportional 
to the classical probability W for an electron with a velocity v and momentum p to return 
to the vicinity of the starting point (more exactly, into the volume of the order v dt (h/p) 2 

important for quantum interference). The probability W is given by the expression^! 

vh 2 r T <t> 

Wcx— J dtP(r,t)\ r=0 , (1) 

where P(r,t) is the probability density. Here we employ the fact that in a diffusive regime 
it is a smooth function of co-ordinates at the scale of mean free path t. To estimate P(0, t) 
we take into account that the electron diffusion is restricted by finite volume of the sample. 
In such a case, we have 



P(0,t)oci$>xp 

" n,nj_ 



n 2 n 2 . 



(2) 



Here D is the diffusion constant, R is the radius of the ring, while d± is its transverse 
dimension. The numbers n, n± have the meaning of quantum numbers for longitudinal and 
transverse diffusive modes, respectively. For a thin ring, d± <C R, only n± = is important. 
One can see that at Dr^/R 2 1 one cannot replace the sum over discrete n, n± in (0) by 
an integral. On the contrary, only n = is important, and W ~ r^A/h. If this quantity is 
small, one can restrict himself with a single return event. 

Now let us concentrate on the case of external magnetic field. In a magnetic field, the 
number n in the expression (||D has to be replaced by n - $/$o, where $ is the magnetic 
flux embedded in the ring, while $o = 7r^c/e.0 It is clear that the quantum contribution is 
maximal if $/$o is close to some integer number no- If the difference n = uq — $/$o = 
we arrive at the same situation as in the absence of magnetic field - only the mode with 
n = no is important. One can expect that this property is also the case at finite |n| <C 1. 
Indeed, for n ^ n 

jsu/ V* A AR2 

5W oc 2^ 



^n D ^/R 2 + 1/7* D 

(for the last estimate we have assumed Dt^/R 2 ^> 1). Consequently, if 5W 1 one can 
neglect the contributions of all the modes with n ^ no to the probability to return. However, 
the corresponding contribution of the mode with n = n Q is not small at Ar^/h > 1. Hence, 
we arrive at the problem of calculation of the localization contribution in the case 

D/R 2 ^>A/n>l/T 4> . 

In this region we can still use a single-mode approximation, but the perturbation theory 
involving a single return event fails. 

The problem in question was addressed by Altland et al.0lil (see also Refs. |I9|-|2T|). The 



authors used the so-called Q-Hamiltonian approach within the framework of the supersym- 
metric method. An intrinsic feature of this method is that one has to cancel out specific 
non-physical contributions. Therefore the supersymmetric approach involves tedious al- 



gebraic calculations. Consequently, the authors ofOiia extensively used computer-algebra 
package. As a result, intermediate equations remain unpublished because, as it was stated, 
the computer printout covers many pages. 

On the other hand, another approach - the so-called replica method - existsS According 
to this approach, one has to replace the system under consideration by N systems identical 
to the original one and at the end tend N — > 0. Usually, after such a procedure one 
arrives at relatively simple expressions. The limiting transition N — > (if done properly) 
automatically cancels out the non-physical contributions that has to be done explicitly within 
the super symmetric approach. 

To take the full advantage of this property one needs a regular procedure to calculate 
the limit iV — > 0. The aim of the present paper is to suggest a procedure of analytical 
continuation of a non-perturbative two-particle Green's function from integer iV to the 
whole complex plane which includes the point N = 0. Such a procedure allows one to 
calculate the limit rather automatically excluding necessity of direct cancelation of non- 
physical contributions. We obtain an analytical non-perturbative expression for persistent 
current in a mesoscopic diffusive ring and compare it with the results ofSEl. 

The paper is organized as follows. In the Section |TJ basic equations for the fluctuation 
of the number of particles, as well as for the persistent current are analyzed. The effective 
action in the single-mode approximation is considered in Section |HI[ In Section |V| the 



particle number auto-correlation function and persistent current are calculated in the non- 
peturbative region. The results are summarized in Discussion. In the following calculations 
we will put H = 1. Then U will be restored in estimates and final results. 

II. BASIC EQUATIONS 

According to0, the main contribution to the persistent current I can be expressed 
through the magnetic flux $ embedded in the ring, as 

rA 8 

I= V^ ((5iV)2) ^' (3) 
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where {(SN) 2 ) is the particle number auto-correlation function, calculated at a given 
value of chemical potential. The latter can be expressed in terms of single-electron Green's 
functions as 



((5N) 2 ) = f° de 1 de 2 K(e 1 ,e 2 ) 



(4) 



where 



*(ei,e 2 ) = ±f d^dr.i^Gf^r^G^r,)) - (QG^n, ri ))(ZGf 2 (r 2 , r 2 )» . (5) 



Here (• • ■) means the usual impurity average. The quantity ([5]) has been calculated in Ref. (LG 
in the limiting case Ar^, <C ft. Our aim is to go beyond this limiting case, namely to calculate 
the correlation function for arbitrary Ar^/h, keeping pp£ ^> h. For this purpose we employ 
the approach by Efetov, Larkin, and Khmelnitskii^l with minor modifications. Namely, we 
will use the so-called Q-Hamiltonian approach within the framework of the replica trick. 
The confined expression for the correlation function K(u) (u = ex — e 2 ) can be written in 



the form (cf. with Ref. g2j) 



K(u) 



N 2 J DQexp(-F) 



dr x dr 2 J DQe- F Tr[AQ{ ri )]Tr[AQ{r 2 )] 



(6) 



where 



dzTr 



D(VQ 



te 



A[Q,A]. 



+ 2 [iuj 



T <t>) 



(7) 



Here A is the vector-potential, [A, B]_ = AB—BA. Having in mind to take into account only 
an elastic scattering by short-range non-magnetic impurities, we can specify Q as 2iV x 2N 
Hermitian matrices, Q 2 = 1, TrQ = 0, iV is the number of replicas, while 

d o\ 

-.1 



A 



V 



/ 



where 1 is the N x N unit matrix. The parameter 1 is introduced phenomenologically. 
We assume that the phase breaking is due to inelastic processes. Followingii, we use the 
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parameterization 

Q = Aexp(W), W 
where B is an arbitrary N x N matrix. 



< o b\ 



-B + 



III. EFFECTIVE ACTION IN SINGLE-MODE APPROXIMATION 



Consider a ring with the radius R and the width dj_ -C R. Consequently, one can take into 
account only the dependence of the matrices B on the angular co-ordinate ip . Expanding 
this dependence into the discrete Fourier series, B = J2 n B n exp(m^), we introduce the mode 
number n. As it has been explained in Section |, only one mode with n = n Q corresponding 
to min(n — $/$o) is important (this assumption will be justified at the end of Section [IV A| ). 
Retaining only this mode and assuming = (1/ ' R)d j 'dip we get WV V W + V V WW = 0. 
Hence, 



1 . d 



V^Q = ^A^-e w = (V^W)W- 1 sinh W = ^ sinh W . 

R o<p R 



Then, sinh 2 W = (1/2) (cosh 2W - 1) can be expanded as a series in 

/ \ 2fc 

W 2k / ,Nfc I \/ B ^0 B n " 



v 



o 



B n B no J 



The item Tr (AQ) can be treated in a similar way As a result, we arrive at the following 
expression for F 



2A 



(8) 



R 2 V" u $o/ ^ v " \ t 4>i 

Now we arrive at an important point. An arbitrary complex N x N matrix B can be 
described by 2 Hermitian N x N matrices. We defined them as 



B = p exp(iip) , B + = exp(—i<p) p. 



(9) 



The quantity F is dependent only on the matrix p. On the other hand, an arbitrary Hermi- 
tian matrix could be diagonalized, the eigenvalues being real. One can immediately observe 



that the integral over p in the expression (^j) with F taken from Eq. (§) for the correlation 
function diverges. This divergence in fact does not occur, because the eigenvalues of p must 
be defined within a finite interval. Indeed, one has to define the variables p in a way to 
obtain one-to-one correspondence between p and Q. On the other hand, one can explicitly 
show that 



Q = Ae 



cos 



B 



B 



- sin VBB+ 
VBB+ 



sinVJiJF 
VBB+ 

cos VB+B 



( 



cos p 



smpe 



■up 



e ltp sin p — cos p T 



\ 



J 



(10) 



where p T denotes the transposed matrix p. Here we have employed the relationship 



p T = ^J e -ip p 2 e ip 



(11) 



which is a consequence of the symmetry properties of the initial replica Hamiltonian (see 
Appendix [A]). It is clear that the matrix Q is a periodic function of p, and one has to specify 
a region at least not larger than 1 period to get one-to-one correspondence. Moreover, to 
get proper analytical properties (damping is the lower semi-plane of the u;-variable) of the 
action F (f|) we have to define the integration limits as (— 7r/2, 7r/2). Finally, the action F 
reads as 



2A 



D 



1 



—tt ( n — — ) Tr sin 2 p + 2 [iuj I Tr cos p 

R 2 V $ / V T <i> 



(12) 



Now let us transform the variables from B,B + to p, u = exp(i<p), the Jacobian being (see 
Appendix 0) 

D(B,B + 



1 „\N 



2(det u^p) 



(13) 



D(p,u) 

We observe that the variables u can be integrated out and canceled with the denominator 
in Eq. (§). 



IV. PARTICLE NUMBER AUTO-CORRELATION FUNCTION 
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A. Eigenvalue representation 



Now we come back to Eq. (Q) . Taking into account that Tr AQ = 2 Tr cos p we see 
that the integrand is dependent only on eigenvalues of p. Hence, we have to transform the 
variables to the eigenvalues and some other ones which could be integrated out both in 
numerator and denominator. This transform is outlined in Appendix y. As a result, 

K _ WJlW^) [Ef = iCos(7rA^/2)] 2 

{IP) SiW em ' 



where {dX} = nj^j 1 dX®\X®\ 2i + N , while 
8(x) = exp 



7lD ( $\ 2 2 /TTX\ 7T . l\ / TlX 

2A^ { n ° - %) 8111 It J " a - e2) -vJ cos It 



(15) 



The expression (|15) contains 3 dimensionless parameters, 



Txfi ir{ei-e 2 ) huDh 2 

7 "AV ^"^A^' £ "^A' (16) 

where n = (n — $/$ )- It is important to keep in mind the following. If max (7, Q) 3> 1, 
only small A are important. Hence, one arrives at the the result, obtained in the framework 
of perturbation theory0'Ei. However, if both 7 and f2 are small one has to sum multi- 
Cooperon contributions that cannot be done in the framework of perturbation theory. There 
is a substantial simplification in the case 

n<l, but £ c =|j^l»l. (17) 

In this case, only one mode with \n Q — $/$o| 1 is important, and it is the case where 
Eq. flTBD is valid. Consequently, we consider the situation where the inequalities (|T7| ) hold, 
but the quantities 7 and Q can be arbitrary. In fact, the mode uq must be considered by 
a non-peturbative way, while the other modes can be treated within the framework of the 
perturbation theory. 
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B. Analytical continuation 



We are not able to calculate the expression (14) for an arbitrary N analytically. Instead, 



we will perform analytical continuation of this expression to arbitrary N, and then calculate 
its limit at N — > 0. 

Let us introduce the quantity 

N-l -i 

Z N = n / dx k xf +N 9(x k ) , (18) 

where 

6(x) = exp [-£ sm 2 (nx/2) - (ifi - 7) cos(7ra/2)] . (19) 
It is convenient to define 

CN = l n Z N = C R + C A : 

where 



3JV-1 



k-— 



JV+l 
2 

3iV-l 



(X 







^ = 1 V In f 1 dxx~ 2k - 1+5 6(x), (20) 
2 iv+i Jo 

5 is a small positive number which later will be put zero. We introduce this parameter to 
keep the important integrals convergent at the limit iV — > 0. The first step is to express the 
sum over k in terms of contour integral over complex k. For this purpose let us take into 
account that the derivative df[27ri(k + 1/2)]/ dk [where f(z) = (e z + 1) _1 ] has second-order 
poles at integer numbers. Consequently, one can express Cr(A) as 



c± \ dk 



1 r k 
F ± = - I dk' In 



1 dxx ±2k '- 1+s 9(x) 



(22) 



The contours C are shown in Fig. |^. These expressions are correct only if other singularities 
except of the poles of / are unimportant. One can show that the function F + (k) has 
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singularities only in the left-hand semi-plane of complex variable k, while the function F~ 
has singularities only in the right-hand semi-plane of k. To prove this statement one has to 
expand the function 8(x) into a Taylor series. For the following, it is convenient to rotate 
the fc-plane by 7r/2 by introducing a new variable k\ = 2irik. The transformed contours C ± 
are shown in Fig. ^. 

Making use of exponential convergence of the integral due to the properties of df/dk, 
we transform the contour integrals to the integrals along the real axis. For simplicity let us 
assume N to be even. As a result, 



Cr(A) 



=F / dk 



( df(k) \ 
I dk ) 



k±mN" 
2m 



k±3iirN\ 



(23) 



Now we are prepared to perform an analytical continuation over N. We have to do it in a 
different way for the functions F for the reason to be discussed later. For this purpose we 
replace iN by ±iVo in the functions F , respectively. Here iVo is a real positive quantity 
which we are going to tend to zero later. Finally, we have 



No 



dk 



'df(k-7iN y 



'df(k-37iN 



2m 



2m 



(24) 



dk ) \ dk 

As a result, the lowest-order term in the Ao-expansion of the function (n is ^ Nq. Finally, 



we get 



C = 2mN( 



dxx- ik/n ~ 1+s 6{x) 



(25) 



The reason of splitting the function into ( R and ( with the replacements — > ±iN 
is as follows. As iV tends to zero, the integration contour comes infinitely close to the cut 
of logarithm functions which enter the expressions for F . Such a situation is not the case 
for any finite N, and it leads to a non-physical pinch which has to be subtracted. Within 
the above mention procedure such a contribution is pure imaginary while the one we are 
interested in is real. The imaginary contributions to F + and F~ have opposite signs. Thus 
the non-physical contribution is automatically canceled out in the sum ( R + ( A . We want 
to note that these terms are of the first order in Nq. They have to vanishes, otherwise the 
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two-particles Green's function would be divergent. In fact, a similar trick has been used 
by Matsubara to formulate the thermal Green's function technique (see e.g., Ref. |2^ ). Let 
us compare our analytical continuation of the function ( to the analytical continuation of 
the two-particle Matsubara Green's function K(Q m ), where Vt m is the external Matsubara 
frequency. In both cases one has to take two functions regular in the upper (retarded) and 
the lower (advanced) semi-plane, respectively. Then one combines the two above-mentioned 
functions into the one having a cut in its complex plane. The physical reason of such a 
splitting into R and A parts is to cancel out non-physical contributions. In the Matsubara 
case the non-physical contributions to K(Q m ) arise in the point Q m — > and cancel out 
after the similar continuation Q m — > iQ of the sum over Q. 



C. Persistent current 



Following^, we express the current according to Eqs. (§),(f|). On the other hand, 



d 2 

K{e u e 2 ) oc e~ c - — — e c 
oei de 2 



Finally, we get 



cA d ( N ( ei = e 2 = 0) 
1 = — —t^- hm -^5 = ~Jony , (26) 

where J = hcD/R 2 & = eD/nR 2 , while 

G= - r dk 92f{k) /pl dx sin 2 ( 7 r 3 :/2) a ;-^- 1 +^(a;, £1 = 0) 

l J-oo dk 2 tidxx- ik /- 1 + s 6(x,n = 0) ' 1 ' 

One can check directly that at 

h HD , . 
A < — < — , h < 1 

T<j, R 2 

the above expressions lead to the expressions obtained in Refs. |T~0| , [T4| . To show that one can 
calculate the integrals with the help of steepest-descent approach to get0'0 

' = -*«. ^ = ^- (») 
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In the region 

7 < 1, £ c > £ > 1 

one can also develop a perturbation theory. Indeed, only small x in the integrals in Eq. 
are important. The physical reason of this fact is the magnetic- field-induced phase-breaking. 
In this region, we arrive at the result 

S = JL^, / = --L£. (28) 
7r z n z nD Tt 6 n n 



This result agrees with the asymptotic result of Ref. [14] for 7 ^> 1, J R? / Dt^ C n <C 1 



Note that the result (p9|) obtained for 7 1 is valid in the region 



R 2 A/hD < n< 1. 

For the case 7 C 1, £ C 1, where the perturbation theory is not applicable, one can put 
9(x, Q — 0) = 1. As a result, we get C/ = 0.21, the current being 

eD 

/ = ~°- 21 ^- (30) 



We observe a maximum at n ~ JR 2 A/hD, the maximal current being 



I m&x ~eJAD/hR 2 . (31) 



Expressions (p9| ) and (p0|) are fully consistent with the curve calculated intllila with a help 
of computer algebraic package. Let us discuss the dependence of the maximal current on 
7 fa h/Ar^. At 7 ^> 1 the perturbation theoryElEl predicts the maximum of the current at 
n ~ JR 2 /Dt < p, the maximal value being 



'—4^. (32) 

Consequently J max oc 7 -1 / 2 at 7 ^> 1, and it 7-independent at 7 <C 1. In this region, the 
persistent current can be estimated also as 



et> F /Ar cl 



max 1/ 4- 

R \ n 
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where r e i is the elastic relaxation time. The quantity Ar e \/h for a typical metal can be 
estimated as (£/ R)(a 2 /A) where i is the mean free path, a is a typical inter-atomic distance, 
while A is the cross-section of the ring. Eq. ([!(]) shows that at Ar^/h ^> 1 the phase-breaking 
time T0 does not enter the expression for the persistent current. Coming back to Eq. ([I]) we 
have to conclude that at t ^ H/A the electronic wave packet does not smear in space. That 
means that the no-mode of the Cooperon is localized in some sense, the localization length 
being of the order JhD/A. Of course, it does not mean complete localization because other 
modes are still under weak localization conditions. 

The range of parameters where the theory above is applicable and leads to non-trivial 
results can be expressed as 

(1/K)(t+/t a ), 

where K ~ (pd±/h) 2 is the number of transverse channels. The left inequality is the criterion 
for a diffusion motion, the first right inequality is just the Thouless criterion S c ^> 1, while 
the last right one is the condition Ar^/h ^> 1. One can see that one needs low temperatures 
to meet the inequality t^/t c \ ^> K ^> 1, as well as samples of very small size. As far as we 
know, no previous experiments satisfy this set of conditions. 

V. DISCUSSION 

As one can see from the preceding sections, the results of the replica procedure being 
complicated for arbitrary integer N are rather simple in the limit iV — > 0. In this limit, the 
non-physical contributions are canceled out automatically, while within the sypersymmetic 
method that has been done this explicitly. An important feature that leads to such a 
simplification is the employed procedure of analytical continuation which has been done 
before direct calculations. Namely, one has to two functions analytical in the upper (lower) 
semi-plane of the complex plane of N, respectively. The proper analytical continuation is a 
combination of these functions. Consequently, in has a cut at = 0. The procedure above 
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allows one to cancel out automatically the non-physical pinch in the two-particles Green's 
function, which otherwise would exist at N = 0. We believe that such a construction 
is important in general for the calculations involving the replica trick. In such a way we 
reproduce analytically and rather simply the results obtained by a computer algebraic 

package. 
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APPENDIX A: EFFECTIVE ACTION DERIVATION 



Here we re-derive the expression (|6|) followingc3 to make clear important symmetry prop- 
erties. Followingii we use the replica trick and introduce a field operators as 

( / + ^ 



* = {V>i,...,Viv,^iV--,^}, * 



Here ipiipf + ipi = 0. The action can be written as 



F = i J (dr) &+(r)(E - H)V + (r) 



E = EI. 



H = [H + U el (r)]I-(-+i5)A 



CO 



(Al) 



Here H is the free-electron Hamiltonian, / is the 2N x 2N unit matrix, while U e \(r) = 
Uq Ylf 1 S(r — Vi), M being the total number of impurities. The first iV rows of Ti. describe 
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an evolution of the retarded Green's functions, while the last N rows describe an evolution 
of the advanced one. The following step is averaging over the positions of impurities. We 
have 



n J -y ex P 

1=1 J v 



M 



r dr 2N , 



M 



(A2) 



dr 2 ' v 

/=i 

Here we have taken into account that only linear terms in and ipf can enter the continual 
integral for the correlation function (Grassman algebra). Because of the same reason, one 
has to allow only for the terms with different / while calculating the product. For a weak 
scattering and in the thermodynamic limit M, V — > oo, , M/V = const, 

£ ~ exp(5/i + iV) , 
MU 



5fx 



V 



dr * + (r)*(r) 



2 

r = | / dr E ^(r)^/(r)^(»-)^(r) . 



(A3) 



Here <5/i is a shift in chemical potential, while g\ = 2MUq/V (go is the coupling constant). In 
a same way as it has been done in Ref. |22|, we introduce an auxiliary scalar field represented 
by Hermitian matrices Q. As a result, the effective ip A interaction can be decoupled as 



exp 



45 



dr J2 ^{r)^ f (r)^(r)^ g {r) 

S+9 



JVQ exp [-Tr J dr ( 




r)Q*(r))~ 


/ T>Q exp 


-Tr fdr^Q*(rj 





(A4) 



Here we have used the definition 2^vg\r e \ = 1. This expression is just the same as Eq. (17) 
fromil. To analyze symmetry properties of the impurity-averaged Hamiltonian let us take 
into account that the initial Hamiltonian possesses the property 7~Lij\ i j <N = 7~t* i+N j +N . Such 
a property has to be kept after impurity averaging and introducing the field Q. In terms of 
Q it reads as iQij = —iQ* +N j +N . Taking into account Eq. fllOD we arrive at the relationship 
(|TT|). The following steps are exactly the same as inH. 
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APPENDIX B: CALCULATION OF THE JACOBIAN 



Let us arrange the columns of the 2N 2 x 2N 2 matrix d(B, B + )/d(p,u) as 



{B\i, . . . , Bni, B12, . . . , Bnn; B n , . . . , B^ N } 



and the rows as 



{Pll, • • • , PNU Pl2, • • • , PN2, ■ ■ ■ , PAWi 
Mil, • • • , UNI] U\2i • ■ • , Un2i ■ ■ ■ , Unn}- 



Taking into account the matrix identities 



dB = dpu + pdu, 
dB + = u~ l dp — ir 1 du u~ 1 p 



we express d(B, B + )/d(p,u) as 



d(B,B+) ( A n A 12 ^ 



d(p,u) 



A-21 A 2 2 



where are N 2 x N 2 matrices. One can show that 



An = u x 1, A 2 2 = -{u V) x u 1 
A 2 i = p ® p ® . . . p , 



A 12 —u®u 



A' 

1 ,0 „.-l 



A' 



(Bl) 



:b2) 



(B3) 



Here x means the Kronecker product, while <g> means direct productil. Making use of the 
identity!! 

det {A x B) = (detAf (det B) q 

(where q,p are the ranges of the matrices A and B, respectively), as well as the Laplace 
expansion of determinant^ we arrive at Eq. 
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APPENDIX C: VARIABLE TRANSFORMATION 



Consider the set of variables which includes eigenstates A( 6 ) and N eigenvectors X^ b \ One 
can see from the definition PikXj® = X^8 ik xf^ that any vector of the type e lx(b) X^ (where 
is an arbitrary phase) satisfies the equation with the same A^ and p ik - Consequently 
one must exclude N extra variables x^ ■ F° r f ms purpose we require the diagonal elements 
X^ to be real. Consequently, the matrix xj^ can be constructed according to the following 
procedure. The first column, x\ l \ contains N — 1 variables xf~\ % ^ 1, while the last (real) 
one, x[^ is calculated from the requirement of normalization. In the next column, x\ 2 \ 
the last N — 2 variables are chosen are independent. The element xf^ is determined by the 
orthogonality of the vectors X^ and X^\ while the last element, X^ is determined by 
the normalization of |X^ 2 )|. The following elements are determined by continuation of this 
procedure. Note, that all the off-diagonal elements are complex ones, so one can consider 
real (Uj^) and imaginary ((V®) parts. In this way we present iV 2 independent elements of 
the matrix p through N eigenvalues A^, and N 2 — N independent variables uj^ and 
From the definition, Pij = £ fc xf h^X^* , one can express p through {U,V,X} as 

^ = X^ k \l + 5 l3 )- l v} k \l-5 lj )], 

-^ = X^ k \l-5 tJ ) + v} k \l + 5 i] )], 
dV> 

d*L = rjWuW + yWyik) + i(v ( k)u ( k) _ (C1) 

Note that the above formulas do not contain summation over repeated superscripts. To 
calculate the Jacobian we arrange the corresponding N 2 x iV 2 transformation matrix in the 
following way. The columns are labeled by iV 2 "old" variables 

{Pil,Pi2, ■ ■ ■ ,PiN}- 

The rows are labeled by N 2 "new" variables 

N-l N-2 1 N-l N-2 1 N 
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Consequently, as follows from Eq. (|C1D , the first (N — 1) rows contain the common factor 
times quantities which depend only on {U, V}. The next (N — 2) lines contain the 
factor \( 2 \ and so on. The last N lines are {A}-independent. As a result, the Jacobian 



can be expressed as lli=i 



A(0 



2(N-i) 



x (some function of {U, V}). This expression has to be 



multiplied by (detp)* = (n£i A w ) , and we arrive at Eq. (|TJ). 
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FIG. 1. The integration contours C± for N 
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FIG 2. The integration contours C± for N 
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